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Abstract 

A set of r non-Hermitian oscillator Hamiltonians in r dimensions is shown 
to be simultaneously diagonalizable. Their spectra is real and the common 
eigenstates are expressed in terms of multiple Charlier polynomials. An al- 
gebraic interpretation of these polynomials is thus achieved and the model is 
used to derive some of their properties. 
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Non-Hermitian Hamiltonians with real spectra are currently being ac- 
tively investigated with respect to their mathematical underpinnings and 
their physical applications. (See the most recent special issues dedicated to 
this topic [H, U Is}.) We present here a set of r non-Hermitian oscillator 
Hamiltonians in r dimensions that all have real eigenvalues. This system 
is also seen to provide an algebraic model for the multiple Charlier polyno- 
mials and is exploited to derive some properties of these special functions. 
The approach is related to the interpretation that was given of the ordinary 
Charlier polynomials in 0, \^ and of the (i-orthogonal Charlier polynomials 
in Q. 

The monic r-multiple Charlier polynomials @, S, 0] Cf(fc) are indexed by 
a multi-index n = (jii, n 2 , • • • , n r ) G N r with length \n\ = n\ + • • • + n r . They 
are orthogonal with respect to r Poisson measures with different positive 
parameters U\, • • • , o T (collectively denoted by a): 



J2cl(k)k l ^ = 0, 1 = 0,1,--- , nj -l (1) 



fc=0 



for 1 < j < r. (In the following we shall omit the suffix a.) They have been 



shown [7|, [l0| to obey the r nearest-neighbor recurrence relations: 

r 

kC H {k) = C H+ ei{k) + (at + \n\)C n (k) + ^2n j a j C H - e - j (k), 

3=1 



kC H {k) = C H+e - r (k) + (a r + \n\)C H (k) +^2n j a j C n ^ e -.(k), 



(2) 



where e} = (0, • • • , 0, 1, 0, • • • ,0) is the j-th standard unit vector with 1 on 
the j-th entry. Subtracting the above relations pair-wise, one finds that the 
polynomials C^ik) satisfy as a consequence 

C H+el (k) - C H+nl (k) + {ai - aj)Cft(k) = (3) 

for all i, j G {1, • • • , r}. 

Let us introduce the Heisenberg-Weyl algebra W(r) associated to har- 
monic oscillators in r-dimensions. H is generated by the annihilation and 
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creation operators and a* (resp.), i — 1, • • • , r, that satisfy the commuta- 
tion relations 

[at, a,j] = [a+, at] = 0, [at, at] = ^ , i, j = 1, • • • , r. (4) 

Denote by \ni, • ■ ■ ,n r ) = \ni) • ■ • \n r ) the normalized simultaneous eigen- 
states of the r number operators iVj = af(ii: 

ajcii \m, ■ ■ ■ ,rii,--- , n r ) = rii | rci, • • • , rii, ■ ■ ■ , n r ) , 

n^GN, i = l,.-.,r, 
(mi, • • • ,m r |ni, • • • ,n r ) = 5 mum ■ ■■S mr<nr . (6) 

Remember that 

[aja, n a,j\ = -6^, [afa, n a^] = S i:j af. (7) 

The algebra W(r) is represented in this number state basis in the standard 
way: 

di \ni, ■ ■ ■ , rii, ■ ■ ■ , n r ) = y/n~ \ni, ■ ■ ■ , rc* - 1, • • • , n r ) , 
\ni, ■ ■ ■ , rii, ■ ■ ■ , n r ) = \/n^+l \m, ■ ■ ■ , rii + 1, • • • ,n r ) . 

Consider now the set of r Hamiltonians Hi, % = 1, • • • , r, defined as fol- 
lows: 

r r 

ifj = ^ a t a 3 + XI a i a t + a i + a i> i = 1; ' ' ' ) r - ( 9 ) 

It is straightforward to see that the multiple Charlier polynomials simultane- 
ously diagonalize the r non-Hermitian oscillator Hamiltonians. To that end, 
form the states 



m = N jr) ^^£L=\n u ...,n r ), keK (10) 
7^ Vni\---n r \ 
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Let us act on \k)) with H^. 



(11) 



n=o Vni\---n r \ 
(\n\ + Oi) \m, ■■■ ,ni,--- ,n r ) 
+ Vrk\ni, ■■■ ,rii - 1, • • • ,n r ) 

r 

^ r) y - 

r 

C H+e .(k) + (<7< + |n|)Q(fc) +^n i (T i C fJ _e;.(fc) 

i=i 

Invoking the recurrence relations (J2D, we thus have indeed 

Hi \k)) = kHi \k)) (12) 

for all i = 1, • • • , r. So although non-Hermitian, the operators Hi have a real 
spectrum given by the non- negative integers, the states \k)) are uniquely 
defined (up to a constant factor) as the joint eigenstates of the operators 
Hi, % = 1, • • • , r with eigenvalues equal to k. Interestingly the Hamiltonians 
H do not commute pairwise. Indeed it is readily found that 

[Hi, Hj] = ai - aj + (<7» - o-j) (13) 
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for all i, j = 1. 



, r. Remark however that 



(a l -a j )\k))=N^Y J 



C H {k) 



V^i! ■■■n r \ 



n=0 

n~i\ni, ■■■ ,rii ,n r ) 



Hence, by property 



/-\ 1 \\ ( 14 ) 

oo 1 

■ (CH+eA k ) ~ Ca+gjik)) \ni, ■ • ■ ,n r ) 

[H i ,H j ]\k))=0. (15) 



The operators Hi, thus on commute "on shell", thereby reconciling the fact 
that they do not commute and yet have common eigenvectors. As such the 
set of Hamiltonians Hi form a "weakly" integrable system. It is also useful to 
observe that the operators Hi can be obtained from the standard harmonic 
oscillator Hamiltonian in r dimensions 

r 

H = ^2 a ^ a j 

by similarity transformations. Using the Baker- Campbell-Hausdorff formula 
e A Ye~ B = A + [A, B] + ±[A, [A, B]] + ±[A, [A[A, B}]] + • • • , (17) 
it is easily shown that 

Hi = SiH S^ (18) 

with 

r 

S, t = e^\[e-^ a K i=l,---,r. (19) 

3=1 

This explains obviously why the spectrum of the Hi are real. 

We shall now show how this framework can be used to derive properties 
of the multiple Charlier polynomials. We shall focus on the step relations 
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they obey and their explicit expression. Suppose we have operators X such 
that 

[H U X\ = -X (20) 

for all i — 1, • • • , r. Then, 

H l (X\k))) = (k-l)X\k)), (21) 

that is, X \k)) is a simultaneous eigenvectors of the Hi with eigenvalue k — 1 
and is hence proportional to \k — 1)). Similarly, if there is an operator Y 
such that 

[Hi,Y]=Y (22) 

for all i — 1, • • • , r, the vector Y \k)) will be proportional to \k + 1)). It is 
immediate to find these operators X and Y from the known lowering and 
raising operators of H , by exploiting the fact that the H^s are related to 
H by the similarity transformations (|T8|) - (|T9|) . Let X = X]i=i a i a *> where 
at are some constants. It is clear that [Hq,Xq] = —Xq. To obtain Xs that 
will obey ( 120 p . one has to determine the constants so that 

S\XqS^ * = S2XQS2 — • • • — S r X()S r (23) 

In this case, it is found that (1231) imposes no constraints on the a,s and hence 
the r operators 

Xj = SidjS' 1 = cij + <jj (24) 

enjoy the property fl20|) . They will be found to imply r step relations for 
the multiple Charlier polynomials. Similarly, using Y = Y^\=i @i a t as the 
generic raising operator of H , the condition analogous to (1241) is found to 
require that all the coefficients $ be equal. Up to a trivial constant factor 
there is thus only one raising operator Y verifying (122"]) and it is given by 



Y = S i iyat]S i =lT i at)+l. (25) 





In order to obtain step relations for the multiple Charlier polynomials from 
the action of Xj, j — 1, • ■ • , r and Y on both sides of (ITOl) . we need to know 
the precise action of these operators on \k)). The key is to relate \k)) to an 
eigenstate of H . 

Let be a state such that 

H \k)* = k\k)\ (26) 
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Take Si \ k)*. Obviously, 

Hi Si \k)* = SiH \k)* = kSi \k)* . (27) 
That is, Si \k)* is an eigenstate of Hi with eigenvalue k. If now, 

S i \k)* = li S 1 \k) m , i = 2,-..,r (28) 

for all i 7^ 1, with 7, some constants, because of (ITHjl . Si will necessarily 
be a common eigenstate of all Hi, % = 1, • • • , r, with eigenvalue fc. Since this 
last property defines |fc)) (up to a constant factor), we could then posit 

\k))=Si\k)*. (29) 

Let us then identify the state \k)) for which (1291) will hold. It is determined 
by the conditions (I26I) and (|28|) . The most general \k)* satisfying (126]) will 
be of the form 

|fc)* = Af< r) Yl d l2 ,.., lr \k-\l\,l 2 ,--- ,l r ), (30) 

'2,'" i^r 

where dt 2i ...j r are some coefficients to be specified from (|25l) and M^ r ^ is the 
normalization factor. The sum in fl30|) is performed over all li G {0, • • • , k} 
such that |/| = (I2 + ■ ■ ■ + £ r ) < k. 

Using the property e A e B = e vA,B ^e B e A which is valid when the operators 
A and B commute with their commutator, one readily finds from (1191) that 

S i = e n - 0i S i e- ai e ai , i = 2,---,r. (31) 

If we choose the coefficients {c^ 2i ...j r } so that 

e- ai e~ a ' \k)* = \k)* (32) 

for all % — 2, • • • , r, condition (|28|) will be satisfied with 7, = e " 1-0 ^. Now 
(|32|) is tantamount to demanding that 

(ai — di) \k)* = 0, i = 2, • • • , r. (33) 

The r — 1 constraints (|33l) are straightforwardly found to imply the following 
r — 1 recurrence relations for the coefficients {di 2i ... t i r }: 



\/lidi 2> ... h +i,...,i T = Vk - \l\d h ,.. ht -,i r (34) 



which are solved by 



V-UUIW...JA'- (35) 



Hence 



W = W£j {k _^J l<-\l>*---,<r), (36) 

~ (r) (r) 

where Mf. = Mf, do,-,o- This normalization factor is readily determined 
from the condition that \k)* be normalized: * (k\k)* = 1. This yields 

M k ir) = 4 (37) 

f 2 

with the help of the multinomial formula. The action of the lowering and 
raising operators of \k)* can now be directly computed. One finds that 

a i \k)*=^\k-l)* (38) 

and also that 

+ 4 + --- + a+)\k)* = y/r(k+l) \k + l}*. (39) 

Although, we have privileged the use of Si in the preceding considerations, 
this choice must not matter and indeed the vector \k)* that has been found 
is seen to symmetrically satisfy the constraints 

(o< - aj) \k}* = (40) 

for all pair i, j — 1, • • • , r. 

We are now ready to complete the derivations of the step relations. Using 
(J23D, (J29D and ([38]), we can now write that 

(a i + o- i )|fc» = |fc-l». (41) 

Using ([TO]) , this yields 

^N^Cnik - 1) = Ni r \c a+el (k) + a 3 C a {k)). (42) 
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To proceed further we need to find . In order to do that, set n = in 
( 142]) and use the recurrence relations (J2J) to observe that 

Cej(fc) = k - ay (43) 



We thus have 



which implies 



with 



N (T) 

= ^ (44) 



^ = -4ttN (45) 



k\r 2 



iVo = (0,---,0|5 1 |0,---,0) = e- CT1 . (46) 

Returning to fl4"2]) . we finally have using (jHJ), the following set of r "back- 
ward" step relations: 

kC H (k-l) = C H+e - j (k) + a J C H (k), j = l,---,r. (47) 

Similarly, for the raising operator, we see from fl25l) . (129]) and (139]) that we 

have 

(af + ■ ■ ■ + a+ + 1) | A;)) = y/r(k + l)\k + 1)) • (48) 

Using again the coherent sum ( ITOl) . with iV^ now given by ( j4"5l) . steps anal- 
ogous to those that led to (147]) bring one to find the "forward" relation that 
(148]) entail: 

n 

C n (k + 1) = C H (k) + ^njC^ik). (49) 

j=i 

The relation (1471) and (149]) can obviously be combined to obtain a difference 
equation for the multiple polynomials C|. 

The Hamiltonian ifo is known to admit a £/(r) invariance algebra, gen- 
erated by the bilinears afaj,i,j = 1, • • • ,r that all commute with H . One 
might be curious to know what is their inference on the discrete integrable 
system on the lattice N r for which the multiple Charlier polynomials provide 
solutions. As for the joint ladder operators X and Y, the common symmetries 
R of the Hi ([Hi, R] — 0, i — 1, • • • , r.) are obtained from the symmetries Rq 
of Hq that have all their conjugates SiRqS^ 1 ,i = 1, • • • , r differ at most by a 
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constant (that does not matter). Resolving those conditions, it is found that 
the common symmetries must be of the form 

Rij = -Rji = Si[(af H h a+)(a; - a^S^ 1 

, , (50) 

= (Cti H hd, + l)(Oi + (Tj - Oj - CTj) 

From (1401) . it is seen that they all annihilate | k)) owing to the factor (a^ — aj) 
in their expression in the first line. Note that Si(ai — a^S^ 1 = (cij + a,; — 
Oj-aj). 

Observe also that these joint constants of motion are in involution: 

[Rij, Rid] = 0- (51) 
By making explicit, the implication of 

Rij \k)) = (52) 
when R^ acts on the r.h.s. of fflUj) . one obtains the relations 

n iCn+e s -e.i (k) ~ rijCfi + ^ s ^g j (k) 

+ (rii - rij +<7i- o-j)Cn(k) + Cft+ei (k) - Cn+ej (k) (53) 

r 

S=l 

which are in fact consequences of (J3]) when combined with the recurrence 
relations (J2]). 

Our final considerations will bear on the explicit expression of the multiple 
Charlier polynomials and on the generating function. Given (I29p . that is 
|fc)) = Si \k)*, and the fact that Si and \k)* are known explicitly, it is quite 
clear that an expression for the multiple Charlier polynomials should follow 
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from working out the action of Si on \k)*. Let us then proceed. 



S 1 = e - ai e- aia i e ai -e- U2a 



2 . . . p u r^r 



Si 



k \l I j ' ' ' j l>r 



a U ai \k- \t\) • JJe _,T W 



(54) 
(55) 



i=2 



E(—a) n i l ' rij\ 
(n--l-V VTT'"^' 

oo 

^V 1 I A; - = ^ |m> (ni| e-^e 01 

ni=0 



(56) 



(57) 



The matrix elements occurring in (157]) can be obtained following techniques 
used in m and |5J . If one sets 



^ n ,k = (k\e-^ e ai \n) 



(58) 



from 



(k\ afate"^ e ai |n> = (k\ e _<Tiai e ai (af a x - a x - triaj" + cti) |n) , (59) 
one readily finds that ip n & obey the recurrence relation 



This implies that 



1 



-e\) n Vn{k)i)Q,k, 



(60) 



(61) 



where p n (k) are polynomials of order n in which are identified with the 
monic Charlier polynomials. The "initial" value ipo,k is directly computed to 
be 



0,fe 



(62) 
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This therefore implies that 



*-|l|)=£ \ 1} = P k -^i) \n x ) . (63) 



m=0 yfm\(k- 



One then uses (155]) . f )56|) and fl63|) to write down the action of Si on \k)* as 
given by fl36|) and fl37|) . First, one calls upon the following expression for the 
standard Charlier polynomials 



Pn(k) = {-aiY 2F0 



-n, — k , 1 
- 01 



^ i! V 



(64) 



where we are using the Pochhammer symbol defined by 

(o) r = o(o + l)---(a + r-l). (65) 

One then identifies the coefficients of |ni,- • ■ ,n r ) on both sides of \k)) = 
Si \ k)* to find that 

c m_ yy fc! (l^|-fc)r , s n 2 \---n r \ 



(-(7- 2 )'»-« (-a,>- s 



(66) 



After some simplification, one finally arrives at the following explicit expres- 



sion: 

ni 



Cn(k) = £ • • • £(-«i)z 1 (-n 2 )i 2 • • • (-n r )i r (-k) h+ ... + i r 



h\ l 2 \ l T \ 

This formula has been given in 0| for the case r = 2. 

Last, a generating function can be obtained from (129]) . by introducing the 
Bargmann representation where the states |ni,n 2 , . . . ,n t ) and the oscillator 
operators cii,af , % = 1, 2, . . . , r are realized by 

K,n 2 , . . . ,w r ) = -= 1 z" 1 z% 2 . . . zl^ , a i = d zv al = z i (68) 

Vni!n 2 ! . . . n r \ 
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Relations OH]) are obviously enforced in this model. 
Recall from Qgg} and (J2HD that 



e?*- n Si\k)*=\k)) (69) 
for any i = 1, 2, . . . , r. In the Bargmann representation 

$ = e -»l e -(ffl*l+-+ffr*r) e a, j (70) 

and 

w = S E 0-$ (71) 

il+«2 + --- + «r=fc 

Since 

) — f( z li z 2i ■ ■ ■ ■> z i + 1? • • • ! 2t)j 

we have, combining ( 16"9~j) and (I7T|) 



"*/2 ^ l x \...l r \ 

h,...l r 1 r 



(72) 



\^ 

n\,...n r 



z^ ... z n J 



This readily gives, the following generating function with the help of the 
multimonomial formula 

e -fa*+™+*T*-\ Zl + z 2 + ■ ■ ■ + z r + l) fc = 

5>(^...^ (73) 

n 

Notice that, as should be, any operator Si, i = 1, . . . , r can be used to derive 
the identity (T73|). 

In summary, let us recall some of our findings. We have presented 
a remarkable set of non-Hermitian harmonic oscillator Hamiltonians in r- 
dimensions with real spectra. Their common eigenfunctions have been seen 
to be given in terms of multiple Charlier polynomials. Had we passed to a 
coordinate representation, we could have obtained relations involving prod- 
ucts of Hermite polynomials and their translates [H, [|5j]. This physical settings 
has provided an algebraic model for the multiple Charlier polynomials that 
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has been used to offer alternate demonstrations of some of their structural 
relations. It would be of interest in our opinion to pursue the algebraic in- 
terpretation of multiple polynomials. The multiple Meixner polynomials in 
particular, should lend themselves also to an oscillator modelization. We 
plan to return to this in a future publication. 
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